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The possibility to avoid the cosmic initial singularity as a consequence of 
nonlinear effects on the Maxwell eletromagnetic theory is discussed. For a 
flat FRW geometry we derive the general nonsingular solution supported by a 
magnetic field plus a cosmic fluid and a nonvanishing vacuum energy density. 
The nonsingular behavior of solutions with a time-dependent A(i(:)-term are also 
examined. As a general result, it is found that the functional dependence of 
A(t) can uniquely be determined only if the magnetic field remains constant. 
All these models are examples of bouncing universes which may exhibit an 
inflationary dynamics driven by the nonlinear corrections of the magnetic field. 

PACS number(s): 98.80.Bp, 11.10.Lm.40.Nr 98.80. Cq 

I. INTRODUCTION 

A fundamental difficulty underlying the standard Friedmann-Robertson- Walker (FRW) cosmology 
is the prediction of an initial singular state where all curvature invariants and some material quanti- 
ties like pressure, energy density and temperature become infinite [1]. Generically, "the break down 
of the laws of physics at a singularity" is a clear manifestation of mathematical inconsistency and 
physical incompleteness of any cosmological model. In this way, although strongly supported by the 
recent observations at low and intermediate redshifts, the present big-bang picture with dark energy 
must be improved at very high redshifts. 
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In order to solve such a problem several attempts based on many disparate mechanisms have been 
proposed in the literature. Some earlier approaches trying to develop a well behaved and more 
complete cosmological description include: quadratic Lagrangians and other alternatives theories 
for the gravitational field [2], a creation-field cosmology [3], a huge vacuum energy density at very 
early times [4], nonminimal couplings [5], nonequilibrium thermodynamic effects [6], and quantum- 
gravitational phenomena closely related to a possible spontaneous birth of the universe [7] . 

More recently, a new interesting mechanism aiming to avoid the cosmic singularity has been 
discussed by De Lorenci et al. [8] through a nonlinear extension of the Maxwell eletromagnetic 
theory. The associated Lagrangian and the resulting electrodynamics can theoretically be justified 
with basis on different arguments. For example, the nonlinear terms can be added to the standard 
Maxwell Lagrangian by imposing the existence of symmetries such as parity conservation, gauge 
invariance, Lorentz invariance, etc [9,10], as well as by the introduction of first order quantum 
corrections to the Maxwell electrodynamics [11,12]. It is worth notice that nonlinear corrections 
may also be important to avoid the black hole singularity. Actually, an exact regular black hole 
solution has been recently obtained with basis on the Einstein-dual nonlinear electrodynamics as 
proposed by Salazar, Garcia and Plebansky [13,14]. Note, however, that the purpose of the present 
work is not to give a detailed description of the possible nonlinear theories and their physical effects. 
Our basic aim is rather modest - we will try to gain new insights into the possibilities beyond Maxwell 
theory based on the simplest nonlinear electromagnetic lagrangian and its connection to the cosmic 
singularity problem. 

In this concern, it has been found [8] that the primordial singularity can be removed because the 
nonlinear corrections reinforce the negative pressure at the early stages of the universe. The authors 
also argued that the nonsingular behaviour is unaffected by the presence of other ultrarelativistic 
components obeying the equation of state p(^ur) — P{ur)/^- 

In the present work, we extend the analysis by De Lorenci et al. [8] using different ingredients. 
We first obtain the behavior of the scale function for a spatially fiat FRW geometry (A — 0), 
thereby showing that their results is a particular case of the general solution. It is also analyzed 
how this solution is modified by the presence of a A-term (vacuum energy density) which may be 
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constant or a time varying quantity. Solutions with constant vacuum energy density also leads to 
a non-singular universe, and the same happens with a time-dependent A(i)-term. However, there 
are singular solutions where the decaying vacuum supplies the energy to a constant cosmological 
magnetic field. In this case, the time dependence of the cosmological term can uniquely be determined 
and corresponds to a slightly modification of the more frequent forms suggested in the literature 
[15-17]. The main physical restrictions, including the time interval where the nonlinear corrections 
must be important are also discussed. 

The article is organized as follows. In section II, we write down the Einstein field equations (EFE) 
for a fiat FRW geometry supported by an energy momentum derived from the extended (nonlinear) 
eletromagnetic Lagrangian. In section III, we generalize the solution derived in [8] which assumes 
a vanishing cosmological term, a spatially fiat geometry and a time-dependent magnetic field. In 
section IV, we obtain a new solution that takes into account the presence of a constant A. The 
behavior involving a time-dependent A is discussed in section V. Finally, in the conclusion section, 
we summarize the basic results and also present some suggestions for future work. In what follows, 
Greek indices run from to 3, Latin indices run from 1 to 3. Unlike of Ref. [8], we adopt the 
International System of Units (see the Appendix on units and dimensions of Ref. [18] for further 
details.) 

II. BASIC EQUATIONS 

As widely known, the Lagrangian density for free fields in the Maxwell electrodynamics may be 
written as 

^(MAXWELL) = —-. F'^Tj^^ = —- F , (1) 

where F^^'^ is the electromagnetic field strength tensor and fio is the magnetic permeability. The 
canonical energy-momentum tensor is then given by 

1 



j:(MAXWELL) _ I J_ 



In this work we consider the extended Lagrangian density to the electromagnetic field 



(2) 
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C = — ^—F + ujF'^ + 77F*^ (3) 
4/^0 

where uj and 77 are arbitrary constants, 

F* = f;,f''- , (4) 

and F*^ is the dual of Ff^^,. As one may check, the corresponding energy-momentum tensor becomes 

T,. = -^^F,-F^^ + (^F* - g,^ . (5) 

Let us now consider the above expressions in the context of a homogeneous and isotropic FRW flat 
model 

ds^ = c'df - a" [t] [dr^ + r'^dd'' + sin^ Q d<\>'\ . (6) 

Naturally, electromagnetic fields may source of the above background only if the fields are considered 
in its average properties [19]. Now, applying the standard spatial averaging process we set 

< F, > = , (7) 

<Si> = 0, (8) 

< E,E, > = —E^Qij , (9) 

< BiB^ > = —B'^Qij , (10) 

< EiBj > = . (11) 

Equations (7) - (11) imply that 

< F,^F% >= I { eoE' + + eoE' - ^] g,^ , (12) 



where t/^ is the four velocity Under such conditions the average value of the energy-momentum 
tensor takes the perfect fluid form, namely: 



<T,,>^{p + p)Mji-pg^^ (13) 



where the density p and pressure p have the well known form 



"4 (-^^^3- (") 

P = 5 P . (15) 



In order to analyze the modifications implied by the use of the modified Lagrangian (3), we 
also assume that for the stochastically defined electromagnetic fields only the average value of the 
squared magnetic field survives at the very early Universe, i. e., we use Eqs. (7) - (11) with 
< E"^ >— Q. In this concern, we remark that a homogeneous electric field in a plasma must give rise 
to an electric current of charged particles and then rapidly decay Indeed, unlike what happens with 
the magnetic field, at present there is no basis whatsoever to presume the existence of an overall 
electric field. Indeed, since the late sixties, it has been recognized that cosmological models with 
an overall electric field bears hardly any relation at all to reahty (see for instance, Zeldovich and 
Novikov [20]). Naturally, this does not means that the < > term appearing in our stochastic 
approach can be neglected in comparison with < B^ >, but one may expect that its influence might 
be small for some special regimes, as for example, when the plasma may be treated with basis on the 
magnetohydrodynamics approximation. For a dense ionized gas, for example, the collision frequency 
can be so high that the electric field and its momenta may arise only as a consequence of the motion 
of the fluid, or as a result of the external charges distribution or (not "frozen in") time- varying 
magnetic fields. Such a behavior may happen in the primeval plasma (below Planck's temperature) 
because the Debye screening radius ~ (T/n)^/^ is very small in comparison with the macroscopic 
relevant scale for nonsingular world models, namely, the Hubble radius cH{t)~^. Keeping these 
remarks in mind, we return to the basic equations by assuming that < E^ > has been neglected so 
that (13) still holds, but the energy density and pressure now read 

P = T^B^ (l - S^io'^B^) , (16) 

P = (1 - AO^iouB') ^\p- ^^uB' . (17) 

D/^o 3 3 

Note that the weak energy condition (p > 0) is obeyed for 



B < ^ } , (18) 
whereas the pressure will reach negative values only if 

B > — , ^ (19) 
On the other hand, there is a widespread belief that the early Universe evolved through some phase 



transitions, thereby yielding a vacuum energy density which at present is at least 118 orders of mag- 
nitude smaller than in the Planck time [21]. Such a discrepancy between theoretical expectation 
(from the modern microscopic theory of particles and gravity) and empirical observations constitutes 
a fundamental problem in the interface uniting astrophysics, particle physics and cosmology, which 
is often called "the cosmological constant problem" [21,22]. This "puzzle" together with the obser- 
vations of type la Supcrnovae [23] suggesting that the cosmic bulk of energy is repulsive and appears 
like a dark energy (probably of primordial origin) stimulated the recent interest for more general 
models containing an extra component and accounting for the present accelerated stage of the Uni- 
verse [24]. A possible class of such cosmologies is provided by phenomenological models driven by 
a constant or a time-dependent A(t)-term (see, for instance, [14-16, 24-30] and references therein). 
The effective time-dependent cosmological term may be regarded as a second fluid component with 
energy density, p^it) = K{t)c'^ /^-kG, which transfers energy continuously to the material medium. 
The conditions under which this kind of cosmology can be described by a scalar field coupled to a 
perfect fluid has also been discussed in the hterature [24-30]. For the sake of generality, we focus 
our attention to a decaying A model, but now in the presence of the primeval magnetic field. In the 
background defined by (6), the EFE read 

while the energy conservation law can be written as 

p + 3-(p + p) = --— . (22) 
a on Lr 

where an overdot means derivative with respect to the cosmic time t. As we shall see, for a constant 
magnetic field it is not need to assume a phenomenological expression for A{t) as usually done since 
it can be uniquely fixed from the above set of equations. 
Replacing (16) and (17) in the Einstein equations (20) - (22) we get 

(1 -8/xoa;52) + — , (23) 



3c^ /lo 3 ' 



— (1-24/^0 , (24) 



a 3 3c^ fjLo 

-(l-16^^ou;B') (b + 2^b)=-^. (25) 

Inserting equations (16) and (20) into (21) we find 

aa + d^-2AeV-f^)5V = 0. (26) 
3 V 3c^ / 



The term proportional to 5^ comes from tlie nonlinear correction in the Lagrangian. Note also 
that if the couphng constant w is zero, the standard FRW differential equation for a radiation filled 
universe plus a cosmological A-term is recovered. By solving any two of the above equations one 
may discuss if the nonlinear terms added to the Maxwell Lagrangian may alter the primeval singular 
state. In order to compare with some previous results presented in the literature, we start our 
analysis by considering some interesting particular cases. 

III. NONSINGULAR MODELS WITH A = 

This is the case studied by De Lorenci et al. [8] where a particular nonsingular solution was 
found. In what follows the corresponding general solution is given and we also show how to recover 
the quoted result. First of all, we remark that if S is a time-dependent quantity and A remains 
constant, equation (25) can be easily integrated to give 

B{t) = B,{^^y, (27) 

where Bq is a constant of integration. In this paper, the subscript does not mean the present day 
value of a quantity. It only indicates its value at an arbitrary time to which appears in the general 
solution of a{t) as a second integration constant. This constant to was arbitrarily chosen in [8]. Thus 
Bo — B{to) if Oo = 0(^0) (in the quoted paper oq was also normalized to unit). We stress that the 
scaling solution (27) holds even for a nonvanishing constant A. Inserting B{a) into (26) one finds 

" I 3^2 1 ^ = 0- (28) 

As one may check, the general solution of the above equation is 



a{t) = ao AaKt- tof + Aaol3o{t- to) + 1 



nl/4 



(29) 



where we have defined 



ao = 



l47iG 



Bo , 



(3o = ±^l-8iiou;Bi . 
In order to compare with the results of [8] we recast (29) in the form 

a{t) ^aoi^alt'^ + Aao^fot + Sq^^^ , 



with 



1o = Po- 2Q;oto , 



(30) 
(31) 



(32) 



(33) 



So = 4q;o^o («o^o - /3o) + 1 ■ 



(34) 



The hnear term in t inside the parenthesis of (32) docs not appear in the solution given by the 
authors of Ref . [8] . This happens because the arbitrary integration constant to was chosen to be 

Po 1 /3/Xoc2(l-8/Xoa;S2) 



^° 2ao 2Bo V ^ttG 
The time behaviour of the magnetic field is readily obtained from (27) and (32). One finds 

Bit) 



(35) 



(36) 



{4alt^ + 4aojot + So)^/^ ' 
with the energy density and pressure defined by Eqs. (16), (17), respectively. Note also that the 
Hubble parameter can be written as 



ao [2ao{t - to) + Po] 



a [4al{t - to)' + 4ao(3o{t - to) + 1] ' 



(37) 



which becomes for t — to 

Ho = H{to)^aoPo^ Bo 



AttG {1-8 Ho ujB^) 
3 fiQ 



(38) 



(In [8], the notation H is used to represent the magnetic field.) From (32) we see that, for large t, 
we recover the classical solution for radiation dominated universes, a{t) oc t^^^. Alternatively, 



we observe from the relation (16) that this solution is recovered for values of t where the magnetic 
field obeys the condition 

SnoOoB'^ « 1. (39) 

The most interesting feature of (32) is that the quadratic function inside the parenthesis does not 
have real roots a; > 0, being positive for any t. Therefore, the model is non-singular with a{t) 
reaching the minimum value amin = 0,0 {8 jj.oUjB^y^'^ at a time tmin = = ^0 — It thus 

follows that the universe is a bouncing one: it begins arbitrarily large at t <^ tmim decreases until 
its minimal value at tmin ^-nd then begins to expand. The values of the magnetic field and energy 
density at tmin are 

p{tmin) = 0. (41) 

Before proceed further, it is worth notice that if (32) describe rightly the evolution of the universe 
in the distant past, it implies the existence of an inflationary era (a > 0) on the interval 

tmin tj <i t <C. tmin ~l~ tl i (^^) 

where 



(43) 

Figure 1 shows the scale factor, the magnetic field, the energy density and the pressure as a 
function of time for a definite value of Bq. From relations (36) and (39) we stress that the classical 
solution is recovered for times much larger than tmin with the universe entering in the standard 
radiation phase. As shown in Fig.l, the nonlinear corrections are relevant only for ^HqujB"^ <; 1/10 

01° t ^ tmin ~l~ 

IV. NONSINGULAR MODELS FOR A ^ 

For constant A, it is easy to see that equation (26) becomes 
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aa + a - -Ac^a^ - I ^ ^ \a =0. 



2 

3^ 



(44) 



V 3c2 

Let us now search for an exact description in the presence of a cosmological constant. By combining 

equations (23) and (27) for a constant A we find 
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where the auxihary scale factor Z and A-term are defined by 



A 



Equation (45) can be easily integrated to give 



Ac2 



1^1/4 p 



Co 



1/4 



where 



(45) 



(46) 

(47) 



(48) 



(49) 



The magnetic field is 

B{t) = 2Bo\f\ 



Do = al (Q;^ + 32A//oa;S^) ■ 



g4VA(t-to) + go g-4VA(t-to) _ 2 ^2 



-1/2 



(50) 



(51) 



It is straightforward to see that the term inside the square brackets of (48) is positive for all t and 
that the scale factor reaches its minimum value 

1/4 



ao 



2A 



+ 32 A /xo ojBq — tto 



(52) 



at 



i'min — ^0 "I 1= ill I 



(53) 



8x/A \Cl} ■ 

As in the previous case, the universe bounces at tmin and, if the solution would effectively hold near 
tjnini an inflationary phase would take place for all values of t such that 
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where 



X = e^^(*-*o) 



The magnetic field at tmin is 

B{t, 



2nG 



+ 



1/2 



^Prom relations (39) and (51) the de Sitter classical solution is recovered for 



t «to + 



In 



and 



t»to + 



1 



In 



(54) 



(55) 



(56) 



(57) 



(58) 



Similarly to what happens with solution (29), the classical solution is recovered for times much 
larger than t„iin- In Figures 2 and 3 wc show the scale factor, the magnetic field, the energy 
density and the pressure as a function of time, for some values of a;, A and Bq. In analogy with 
the solution (32), we have that the nonhnear corrections are relevant only for S/jLoloB^ > 1/10 or 

21 



t ^ to 



+ 



In 



1 + A 1- 



V. NONSINGULAR MODELS FOR A TIME-DEPENDENT A 



The possible cosmological consequences of a decaying vacuum energy density, or A(^) cosmologies 
are still under debate in the recent literature [14-30]. Such models may also be described in terms 
of a scalar field coupled to a fiuid component. Another important motivation is its connection with 
the cosmological constant problem. In general grounds, one may expect that a decaying vacuum 
energy must play an important role on the universe evolution (mainly in the very early Universe) 
and, probably, more interesting, it may indicate suggestive ways to solve the A-problem, as for 
instance, by describing the effective regimes that should be provided by fundamental physics. In 
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the majority of the papers deahng with a time- varying A, the decaying law is defined a priori, i.e., 
in a phenomenological way. The most commonly postulated decaying laws are those in which h.{t) 
decreases as some power either of the scale factor a{t) or the Hubble parameter H (see [17] for a quick 
review) . Some authors have also considered scaling laws formed by a combination of both quantities 
[16]. As remarked before, these proposals are in accordance but do not explain the difference of more 
than 100 orders of magnitude between the cosmological constant value at the beginning of universe 
(provided by particles physics) and its actual value estimated from cosmology. In general, the EFE 
imply that once h.{t) is given one may integrate them for obtaining B{t) and a{t). Conversely, for a 
given dependence of -B(t), a unique time dependence for A(t) is readily fixed by the field equations. 

Let us first analyze phenomenological models with a cosmological term defined by A = 3/3c~^i?^, 
where /3 is a positive parameter smaller than unity [16,30]. The differential equation driving the 
scale factor is readily derived by combining relations (16), (17), (20) and (21). One finds 



ij + 4(1 -/?)//'- (1 - /?) 



1 - 



4(1 -/9^) 



an 



(59) 



where H — a/a. The constants ckq is defined by relation (30) whereas /3* is given by 
Now, separating the variables in (59) one finds 



(60) 



an 



A simple integration results 



1 - 



1 4(1-/3.^) ^, 



al 



-AH^\ dH = {I- (3){t-to). 



H 



a 
a 



[4«^(1 - py{t - toy + 4ao/?*(l -P){t- to) + 1] 
and integrating again, we obtain for the scale factor 

a{t) = ao kalil- P)\t-tof + 4 ao{l - P) P^t - to) + f f'^^^'"^^' 



(61) 



(62) 



(63) 



This solution is nonsingular for w > 0, with a{t) reaching its minimal value, a^m = (8a*o<^(1 ~ 



^)^2)i/4(i-/3)^ at a time t. 



tn - 



It thus follows that the universe is a bouncing one. It 



2ao(l-/3)- 

begins arbitrarily large at i -C ^mm; decreases until its minimal value at t^in and then begins the 
expansion phase. For completeness, the expression for the magnetic field is 



12 



Bit) = ^o(^ - f^y^' (64) 

and at tmin, it is readily checked that B{tmin) — 2^2^Mot^ p{^min) — 0. As one should expect, in 
the limit = 0, all the results above for a time dependent A-term reduce to that ones of A = 
(see equations (27)- (34)). Before proceed further, it is worth notice the existence of an inflationary 
era (a > 0) which depends on the value of the (3 parameter. For 1/2 < /3 < 1 the universe always 
evolves trough an accelerated expansion state. However, if < /5 < 1/2, it inflates on the interval 

tmin -ti < t < tmin + tj, whcrC 



tl^ 



\ 71(7(1-2/3) ■ 

In Figures 4 and 5 we show the time dependence of the scale factor, A-term, magnetic fleld, 
energy densities and pressure for some selected values of Bq and /3. The main conclusion is that the 
singularity must be avoided for a generic time dependent A. As shown in Figs. 4 and 5, the nonlinear 



corrections are relevant only for ShqluB^ > 1/10 or t < tmin + 3y8/Xo'^(1 ~ /5)-So/2q;o(1 ~ P)- 

At this point one may ask by the inverse treatment, i.e., if the singularity is avoided for a given 
B{t). Such a question is immediately answered by examining the simplest case, namely, that one for 
which the magnetic field remains constant in the course of the evolution. This possibility is clearly 
allowed by the generalized energy conservation law (see (25)). The energy density of the magnetic 
field is kept constant because energy is continuously drained from the decaying vacuum component 
to the magnetic field. Actually, if A-term is maintained constant the unique solution with a constant 
magnetic field is the trivial one (B = 0). Unlike the previous solutions, we note that there is no an 
analogous classical solution for constant magnetic field. Therefore, we will analyze such a possibihty 
regardless of any constraint on its domain of validity. 
If B{t) = Bq = constant, the energy conservation law yields 

A{t) = Ao + ^Koln(-) , (66) 

where Aq = A (to) and 

IGttG B^ , , , 

Ko = (1 - 16/xoc^5o') ■ (67) 
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Substituting (66) into (23), we get for the scale factor 



a{t) = Oo exp 



= Co exp 



it - toY + Ho{t- to) 



(68) 



where 



\ 3 3c^ iio 



Pi = - Io-Hq] and /3o = —r- to - Ho 



(69) 



(70) 



The Hubble parameter is 



and we have H = for 



H{t)^^{t-to) + Ho^^t + P^, 



2P, _ 2Ho 

tc — — TT — o — ^0 — 



(71) 



Koc^ Koc^ 
At this point, the scale factor reaches the value 



(72) 



(73) 



The behaviour of the solution will depend on the sign of the constant Kq (for Ko = we get the 
de Sitter solution). For Kq > , that is, for Bq > 1/(4^//!^), the universe is always accelerated 
(a > 0) e has a minimum size at tc- 

A much more interesting solution is the one corresponding to /Tq < (-Bq < l/(4^/Xo<^))- For this 



range of Bq, a{t), a > ior t < tc — ^J—2Jk^ and t > tc-\- ^j—ljKoC^ and a{t) has a maximum 
at tc- It is worth notice that the time interval At(7v/), prior to tc-, for which the solution is not 
inflationary depends on the value of Bo as 



At (AT/) = 2< 



Koc^ Bo\^i:G{\-\^iio^Bl)' 
The cosmological term A dominates the dynamics of universe for values of t where 

Ac^ > STrGp. 



(74) 



(75) 
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For both models {Kq > and Kq <0), the condition (75) is satisfied hy t < and t > where 



and 



-2Hoc + 2^/^f^ 



^4 = r.J ' ■ (77) 



In Figure 6, we show the scale factor and the cosmological term as a function of time for Kq > 
and some values of Bq and Aq = Aoc^/3. The same quantities have also been plotted for /Tq < in 
Figure 7. 

Naturally, if one expects any such model to properly describe the evolution of the real universe, it 
would be advisable to take into account other matter fields, such as ultrarelativistic matter, scalar 
fields or dust. In [8], it was demonstrated, for the case A — 0, B — B{t), that the presence of ultrarel- 
ativistic matter with an equation of state P[ur) = P(ur)/3 would just amount for a reparametrization 
of the constants Bq and lu. 

At this point, we would like to stress the mathematical consistence of the whole set of solutions 
derived in the present work. In general, there are 3 unknown functions: the scale factor a, the 
magnetic field B and the cosmological term A (constant or time-dependent). As one may check for 
each case, the number of unknown functions and equations coincide, the unique exception is related 
to models containing a variable A{t) term for which the phenomenological law A = 3/3c~^i?^, has 
been considered (see Refs. [16,18]). 

VI. CONCLUSION 

We have examined whether nonlinear corrections to the Maxwell electrodynamics may avoid the 
cosmic singularity occurring in fiat FRW universes. In brief, the answer is positive. We show that by 
discussing a large class of analytical cosmological models under three different assumptions. In the 
first case, the cosmological A-term is identically zero and the dynamics is driven by a time dependent 
magnetic field. This class generalizes the particular solution previously found by De Lorenci et al. [8], 
and confirms their statement concerning the avoidance of the initial singularity. In principle, since 
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the solutions are non-singular, they potentially solve the horizon problem. We have also examined if 
the basic features of such models remain true if new ingredients are introduced in the matter content. 
In this concern, models with a constant and time dependent A-term were studied with some detail. 
Again, for both cases, the universe is also non-singular, bouncing at a critical time when the scale 
factor reaches its minimum value. 

For a decaying vacuum energy density we discuss two different scenarios. In the first one, it was 
phenomenologically described by A{t) ~ as assumed by several authors [16,17]. These models 
are nonsingular and resemble the solutions with no A. The second scenario is a rather curious 
solution which describes a universe driven by a constant magnetic field. The time behaviour of the 
cosmological term is now uniquely determined by the EFE as a logarithm of the scale factor. It 
should be interesting to examine if such results are maintained in the presence of other matter fields, 
as well as for universes with non-zero curvature. 

Finally, in analogy with the cosmological case, one may ask if nonlinear terms in the Maxwell 
Lagrangian may remove the physical singularity present in a charged black hole (Reissner-Nordstrom 
solution). This problem will be discussed in a forthcoming communication. 
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a„ t 



FIG. 1. The upper panel shows the scale factor (solid line), the magnetic field (dashed line) and the classical solution (doted 
line) for (w = 0). The lower panel shows the energy density (solid line) and the pressure (dashed line) for the model with A = 
and 2Bo^fEjii^ = 0.2. 




a„t 



FIG. 2. The upper panel shows the scale factor (solid line) and the magnetic field (dashed line). The lower panel shows the 
energy density (solid line) and the pressure (dashed line) for the model with a constant non-vanishing A . The values for A 
and Bo are such that y/X/ao = 0.4 and 2Bo^/2iJJIo = 0.2. 
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0.0 



FIG. 3. As in Figure 2 but for VT/ao = 0.01 and 2BoV2a^ = 0.1. 




FIG. 4. The upper panel shows the scale factor (solid line), the magnetic field (dashed line) and the classical solution (doted 
line) for (a; = 0). The lower panel shows the energy density of magnetic field (solid line), the energy density of the A term 
(dashed line) and the pressure (doted line) for the model with A = ^H^, P = 0.4 and 2So\/2cJ/Uo = 0.2. 
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a„(l-P) t 



FIG. 5. The upper panel shows the scale factor (solid line), the magnetic field (dashed line) and the classical solution (doted 
line) for (a; = 0). The lower panel shows the energy density of magnetic field (solid line), the energy density of the A term 
(dashed line) and the pressure (doted line) for the model with A = ^H^, f3 = 0.6 and 2Bo\/2cJ/Uo = 0.2. 
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FIG. 6. The scale factor (solid lino) and the cosmological term (dashed line) for the model with constant magnetic field, 
time-dependent A, i^o > {2Bo\/2iJJIo = 1)- In the upper panel \/Ao/ao = 1 and the lower panel is for s/Xo/ao = 0.5. 
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